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HEAT INVARIANTS OF THE STEKLOV PROBLEM 

lOSIF POLTEROVICH AND DAVID A. SHER 

Abstract. We study the heat trace asymptotics associated with the Steklov eigenvalue 

. . . problem on a Riemannian manifold with boundary. In particular, we describe the struc- 

CO I turc of the Steklov heat invariants and compute the first few of them explicitly in terms 

of the scalar and mean curvatures. This is done by applying the Seeley calculus to the 

Dirichlet-to-Neumann operator, whose spectrum coincides with the Steklov eigenvalues. As 

an application, it is proved that a three-dimensional ball is uniquely defined by its Steklov 

O ' spectrum among all Euclidean domains with smooth connected boundary. 

< 

1. Introduction and main results 

Ph I 1.1. Steklov eigenvalue problem. Let fi be a smooth compact Riemannian manifold of 

^ ' dimension n with smooth boundary M = dVl of dimension n — 1. Consider the Steklov 

'^ ■ eigenvalue problem on fi: 

Am = in f2, 



c^ 



. 1.1.1 

^ = au on M. ^ ^ 

The spectrum of the Steklov problem is discrete and is given by a sequence of eigenvalues 
> ■ 

rn ■ = (To < 0"! < (T2 < ■ ■ ■ /^ OO. 

^ ■ The restrictions of the corresponding eigenf unctions to the boundary form an orthogonal 

l> ■ basis in L'^{M). 

tJ" . Geometric properties of Steklov eigenvalues on Riemannian manifolds have been actively 

^ I investigated in the recent years, see |CEG[ IFSlj IFS2[ IGP2[ iJal IKKPj . In particular, various 

estimates on eigenvalues and their multiplicities have been obtained. 

1.2. Steklov heat invariants. The present paper focuses on geometric invariants associ- 
^\ . ated with the Steklov problem. As follows from Weyl's law for the distribution of eigen- 

j^ I values, the (n — l)-dimensional volume of the boundary M can be determined from the 

Steklov spectrum. In order to obtain further geometric information, we consider the heat 
trace asymptotics. While such an approach is quite standard in spectral geometry (see [GI] 
and references therein), to our knowledge, it has not previously been systematically applied 
in the context of the Steklov problem. 

Steklov eigenvalues an niay be also viewed as the eigenvalues of the Dirichlet-to-Neumann 
operator V : C^{M) ^ C^Im), given by Vf = d^{Hf), where / G C°°{M), Hf e C'^(fi) 
is its harmonic extension to fl and d,^ denotes the outward normal derivative. It is well 
known that D is a pseudodifferential operator of order one (see |Tay[ pp. 37-38]). As follows 
from the results of |DGt |Ag[ IGrSej . the trace of the associated heat kernel, e~*^, admits an 



asymptotic expansion 



oo „ oo oo 

^g-t<x. _ Tre-*^ = / e-''^{x,x) dx ~ ^afcr"+i+^' + ^fe^t'logt. (1.2.1) 

i=0 •^'^^ k=0 1=1 
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The coefficients a^ and 6/ are called the Steklov heat invariants, and it follows from fll.2.ip 
that they are determined by the Steklov spectrum. The invariants Oq, . . • ydn-i, as well as 
bi for all /, are local, in the sense that they are integrals over M of corresponding functions 
ak{x) and bi{x) which may be computed directly from the symbol of T. However, a^ is not 
local for k > n [Gl ]IGiGr] . One of the main goals of this paper is to investigate the local heat 
invariants Qq^x), . . ., a„_i(x). In particular, we discuss their general form and also compute 
Oo, tti, and a2 explicitly. 

1.3. Structure of the heat invariants. We will prove that the Steklov heat invariants 
ak{x) must consist of explicit polynomials in the metric, its inverse, and their derivatives 
at X. For any term in such a polynomial, we let the (total) weight be the total number of 
derivatives in that term and let the normal weight be the total number of normal derivatives 
in that term; for example, a term of the form g°'^'^gap,n has weight 2 and normal weight 1. 
Note that by a Taylor series argument (see Section 2.1), all derivatives of the inverse metric 
at X may be expressed as polynomials in the derivatives of the metric itself at x, with the 
same weight and the same normal weight as the original term. Moreover, we say that a 
polynomial in the metric, its inverse, and their derivatives has total weight k if each term 
has total weight k. We also say that such a polynomial has normal parity 1 if each term has 
odd normal weight, and normal parity if each term has even normal weight. The following 
result is proved in section 13.11 

Theorem 1.3.1. For each k with Q <k < n — 1, the pointwise Steklov heat invariant a^lx) 
is a universal polynomial in the metric and its derivatives, of total weight k and normal 
parity equal to k mod 2. 

1.4. Explicit expressions. We now give expressions for the first three Steklov heat invari- 
ants. Let Hi be the mean curvature, and let H2 be the second order mean curvature, given 
in terms of the principal curvatures Ai, . . ., A„ by 

H2 = 7 TTT rr 2, -^aA/3. 

In — In — 2 -^^ 

Further, let i?^ be the scalar curvature of the domain fi and Ru be the scalar curvature of 
the boundary M. Finally, let 

K, = Vol(§"-^) = ^. 

Theorem 1.4.1. The first three pointwise heat invariants of the Steklov spectrum are given 
by 

ao(x) = (27r)-"+V„r(n-l); (1.4.2) 

ai{x) = ^^^^^ H,, (1.4.3) 



a2[x 



VnTin - 2) 



K27r) 



n-l 



{n - l){n - 2){n^ - n - A) o n{n - 3){n - 2) 

ill H' 

n+1 ^ n+1 



2 



n — 2 n — 4 

H Rq — — ; zRm 

n-l 3(n-l) 



1.4.4) 
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These formulas for afc(x) are valid whenever k < n — 1. 

Remark 1.4.5. Edward and Wu have shown in |EWj that ai = for any simply connected 
domain fi in M^; this agrees with our resuh, which indicates that in fact the pointwise heat 
invariant ai{x) = whenever n = 2. 

Theorem 11.4. II is proved in sections l4.ll and [ 4.2[ The expression for a2{x) simphfies in the 
special case where the ambient space Q has constant sectional curvature K. In this case, 
R^ = n{n - 1)K and Rm = {n - l)(n ~ 2){K + H2) (see [ALMj l We therefore have 

Corollary 1.4.6. When Q has constant sectional curvature K, 



a2{x) 



V„T{n 



K2< 



n-l 



(r2-l)(n-2)(n2-r2-4) 2 4(n2 - 3n - 1) 

n -"1 ?r^ — r\ — -^^^ 

n + 1 3(ra^ — 1) 

2n{n- l){n-2) 



n + 1 



K 



:i.4.7) 



In particular, when n = 3 and f2 C M (i. e. K = Q), 



1.5. Applications to spectral rigidity. The coefficients a^ in the heat trace expansion 
f ll.2.ip are given by 



Ofc = / ak{x} 



As was mentioned in section 11.21 it follows from Weyl's law that the volume of M is deter- 
mined by the Steklov spectrum. Formula (11. 4. 21) provides an alternative proof of this result, 
as 

ao = ^^!_r^^ , Vol(M). (1.5.1) 

2«-%Vr(^) 

For n = 2, the first heat invariant ai = 0, but when n > 3, ai is a nonzero multiple of the 
mean curvature. As a consequence, we obtain 

Proposition 1.5.2. When n > 3, the total mean curvature j^^ Hi of the boundary M = dQ 
is determined by the Steklov spectrum. 

In sections 15.11 and 15.21 we prove the following spectral rigidity result: 

Theorem 1.5.3. Let Q G M.^ be a domain with connected and smooth boundary M. Suppose 
its Steklov spectrum is equal to that of Bp, a ball of radius p. Then Q = B^. 

The proof involves two components. First we use arguments related to the Weyl asymp- 
totics to show that M must be simply connected; this step is based on an argument of 
Zelditch [Z]. We then use the heat invariants we have computed to show that in fact M 
must be a sphere. 

Remark 1.5.4. An analogue of Theorem 1 1.5. 31 in two dimensions is well-known. Indeed, it was 
shown by Weinstock ( [Wej . see also |GP1] ) that among simply connected planar domains, of 
a given perimeter, the first Steklov eigenvalue attains its unique maximum on a disk. Since 
perimeter is spectrally determined, this implies that a disk is uniquely determined by its 
Steklov spectrum among all simply connected planar domains. An alternative proof of this 
result could be found in [EdJ . 
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It is unknown whether a ball maximizes the first Steklov eigenvalue in higher dimensions 
among all domains of fixed volume of the boundary; therefore, the above argument can not 
be used in higher dimensions. There exists a generalization of Weinstock's inequality due to 
[Brj . but it involves both Vol(M) and Vol(f2), and it is not known whether the latter is an 
invariant of the Steklov spectrum. 

We conclude this section with the following conjecture, motivated by Theorem 11.5.31 and 
Remark 11.5.41 

Conjecture 1.5.5. A ball in M" is uniquely determined by its Steklov spectrum among all 
n-dimensional Euclidean domains. 

Note that the analogue of this conjecture is known to be true for Dirichlet and Neumann 
eigenvalue problems. 

There are other situations in which similar rigidity theorems can be proved using heat 
invariants. For example, the round sphere is known to be determined by the spectrum of 
the Laplace-Beltrami operator in dimensions ra < 6 |Tan] . Hassell and Zworski also use heat 
invariants to prove resonant rigidity for 5*^ in the setting of obstacle scattering in M? |HZ] . 

1.6. Plan of the paper. In sections |2. II and [2^2] we compute the symbol of the Dirichlet-to- 
Neumann operator in boundary normal coordinates, following |LU] . Combining these results 
with the Seeley calculus, we obtain local expressions for the first n — 1 heat invariants in 
section 12.31 This allows us to prove Theorem 11.3.11 in section 13.11 Some applications and 
extensions of this result are presented in section 13.21 In particular, the heat invariants of 
product-type manifolds are described in Corollary 13.2.21 Sections 14.11 and 14.21 are devoted 
to the proof of Theorem 11.4.11 In section 14.31 we verify these results by a direct calculation 
of the first few heat invariants of Euclidean balls in dimensions three and four. Finding an 
explicit expression for the heat invariant 02(3;) is the most technically involved part of the 
paper. The proofs of auxiliary lemmas used in this calculation are presented in sections 16.11 
and 16.21 Some integrals appearing in the computations are given in section 16.31 Finally, in 
sections 15.11 and 15. 2[ we prove Theorem 11.5.31 

Acknowledgments. The authors would like to thank P. Hislop and P. Perry for useful 
discussions. Research of LP. was supported in part by NSERC, FQRNT and the Canada 
Research Chairs Program. Research of D.S. was supported in part by the CRM-ISM post- 
doctoral fellowship. 

2. Seeley calculus for the Dirichlet-to-Neumann operator 

2.1. Specialized coordinates. To compute the first few terms in the symbol of P, we 
follow the work of Lee and Uhlmann jLUj . Since the heat invariants ak{x) are local, we are 
free to choose coordinates which are well-adapted to our analysis near a fixed point P and 
then compute ak{P)- We use Riemannian normal coordinates centered at P, with any choice 
of initial orthonormal frame, to write the metric on M in coordinates as 

n-l 

y gap{x)dx°'dx^ , 

o,/3=l 

where P corresponds to the origin and Qapix) = Sap + O(x^). Let the inverse metric be g°'^. 
Notice that all first derivatives of the metric and its inverse vanish at P. 
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Following |LU] . we now extend these coordinates to 'boundary normal coordinates' in 
a neighborhood oi M G fl. For each x & M near P, let x"" be the parameter along the 
geodesic starting at x with initial direction given by the inward-pointing boundary normal. 
Then the coordinates (x"|^z]^, x"), which we write x' = (x, x") (note that the notation differs 
slightly from that in |LUj ). are good coordinates on f] in a neighborhood of P. Moreover, 
the definition of x" does not depend on the initial orthonormal frame. In these coordinates, 
the metric is precisely 

n-l 

g= J2 9ap{x')dx''dx^ + (rfx")l 

a,l3=l 

As in |LU] . we use the Greek indices for coordinates along M and Roman indices for co- 
ordinates in VL. Let the dual coordinates to x" be ^", and let the squared volume element 
5{x') = det{gij{x')) = det{gafi{x')). 

Using these coordinates, we still have freedom to choose our initial orthonormal frame on 
M; in particular, we may choose a frame that diagonalizes the second fundamental form. 
First fix any orthonormal frame Xi, . . . ,Xn-i, and let the local normal be X„ (which is 
independent of the choice of frame). Then we may write the linear operator associated to 
the second fundamental form as 

n-l 

s{x^) = {-Vx.Xrf = 5^-r°„(P)x,. 

a=l 

Using the formula for the Christoffel symbols (see, e.g. |dC] ) together with the fact that 
g{P) = Id, we have 

1 

2' 
However, since gna = in a neighborhood of P in Q, we have 

1 

So the (a, /3) entry of the matrix of the second fundamental form at P is just —^gai3,n{P)- 
This is a symmetric matrix and hence orthogonally diagonalizable, so we may choose an 
initial orthonormal frame which diagonalizes the second fundamental form at P. In these 
coordinates, gap,n{,P) = whenever a ^ (3, and gaa,n{P) = — 2Aq, where Aq is the corre- 
sponding principal curvature at P. 

It will also prove helpful to have similar expressions for the inverse metric at P. Viewing 
the metric g{x') as a function of x" with values in symmetric 2-tensors on M, we Taylor 
expand around x" = 0: 

g{x') = Id +Aix" + A2(x")2 + C((x")3). 

The inverse metric then must have the expansion 

g~\x') = Id-Aix" + {Al - A2)(x")2 + 0((x")'). 

We conclude that (7°^^'" = —ga^^n, and in particular that ^°/^'"(P) = 2Aq(5q^. Moreover, from 
the second order terms, we observe that 

^^«.,nn(p) = ^8A^ - 5^^«.,„.(P). (2.1.1) 



■r^ri — —■^{gna,l3 + gal3,n — fi'/3n,a) 
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2.2. The symbol of the Dirichlet-Neumann operator. Letting D^j = —id^j, we write 
the Laplacian on Q as 

Dl„+tE{x')D,n+Q{x',D,), 

where 

E{x') = -^J29"^^^'^^--9ap{x'), 

a,l3 

Q{x,D,,) = J29"^i^')D.cD,,,-zJ2il9"^i^')d^-'^og6{x') + d,cg»^{x'))D,,. 

a,l3 a,l3 

In particular, we can read off the full symbol of Q, which we write q2{x', ^) + qi{x', ^). 

By observing that the Dirichlet-to-Neumann operator V satisfies a Riccati-type equation, 
Lee and Uhlmann compute its full symbol |LU] . First, define the symbol of a pseudodiffer- 
ential operator with parameter Xn, which we call V, recursively by the following formulas: 

(2.2.1) 



'^o 



— f^rf^v^D 






(2.2.2) 



and for m > 0, 



' — m— 1 



/ 



2v^ 



E 



— m<j<l 
-m<fc<l 

\\K\=j+k+m 



K\ 



df{f-j)D^,{f-k) + 9,nf_„ - Ef. 



\ 



J 



(2.2.3) 



Let T) be the restriction of P to a;„ = 0. Then as in [LUJ the Dirichlet-to-Neumann operator 
T) is equal, modulo smoothing operators, to — P; the sign is chosen so that T> has positive 
principal symbol. Write the symbol of P as ri + tq -|- r_i + 6„2, where 6_2 has order —2; 
expressions for ri, tq, and r_i may be computed from (I2.2.ip . (I2.2.2p . and (I2.2.3p . 

2.3. Seeley calculus and the heat invariants. One may now apply the work of Seeley 
j5e] to compute the local expressions for the first n — 1 heat invariants of V. We follow the 
exposition of Gilkey and Grubb |GiGrj . Namely, let S{X) be a parametrix for V — \: that 
is, a pseudodifferential operator of order —1 with parameter A for which {T> — X)S{X) and 
S{\){V — A) are both pseudodifferential operators of order — oo. Such a pseudodifferential 
operator must have symbol s_i(x, ^, A) + s_2(a;, ^, A) + . . . given by: 

s_i(a;,^,A) 



ri — A' 



/ 



S-l- 



(x,e,A) 



ri — A 



\ 



E 



\\K\ 



— »n<j< — 1 
-m<fc<l 

\|_ft:|=m+fe+j>0 



K\ 



-dfnd^s, 



(2.3.1) 



(2.3.2) 



For later reference, we write out the formulas for s_9 and s^ 



S-2 = -roin - A) - i{d^ri ■ d^ri){ri - A) 



-3. 



(2.3.3) 
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s_3 = -(ri-A)-i 



^0-5-2 + r-is_i - i(%i ■ d^s_2 - di-ro ■ 9^.s_i) - ^ —d^rid^s_i 



K\ 

\K\=2 



(2.3.4) 
The heat invariants may now be computed by using the functional calculus; see [Se] for the 
details. Letting F be a contour around the positive real axis, we have that for < A; < ra — 1, 



ak[x) 



(27r: 



-/ [e-^s,,.k{x,^,X)dXd^. (2.3.5) 



3. General form of the heat invariants 

3.1. Structure theorem. Here we prove Theorem II. 3. 1[ The proof proceeds by first ana- 
lyzing the inductive formulas for the symbol of V and then passing to the inductive formulas 
for the heat invariants. In these formulas, we repeatedly encounter expressions of the form 

\^\-'p{x',0, (3.1.1) 

where fc G Z, |^p is shorthand for g'^^^'^^^-, and p{x',C,) is a polynomial in ^ with coeffi- 
cients equal to polynomials in the metric, its inverse, and their derivatives. We say that an 
expression of the form (I3.1.ip has total weight w if each term in p{x',^) has total weight 
w. Additionally, we say that such an expression has total parity equal to 1 if each term in 
p{x', C,) has odd total weight, or equal to if each term has even total weight; note that the 
total parity is not defined if there are terms of both odd and even total weight in p{x',^). 
Analogous definitions may be formulated for normal weight and normal parity. Finally, we 
define the ^-parity of (I3.1.ip to be 1 if p{x', ^) is odd in ^ and if p{x', ^) is even in ^. 

Note that each of these definitions is independent of a change of form in (13.1.11) - for 
instance, increasing A; by 2 and multiplying p{x',^) by |^p changes nothing. Furthermore: 
suppose p{x',^) = \^\~^p{x',^) has well-defined total weight, normal parity, and ^-parity. 
Then it is an easy calculation to see that each of d^^p, d^^p, and d^ap may also be written 
in the form (I3.1.ip . and each has well-defined total weight, normal parity, and .^-parity. All 
x'- derivatives increase the total weight by 1, the a;"-derivative also fiips the normal parity, 
and any ^-derivative fiips the ^-parity. 

Lemma 3.1.2. For each m > —1, f_m{x') is a sum of terms of the form Ii3.1.1\) . each of 
which has total weight m + 1. Moreover, the total parity, normal parity, and ^-parity of each 
term are all well-defined, and their sum is equal to zero mod 2. 

Note that since V = —V\xn=o, an identical statement obviously holds for rj{x,^). 

Proof. The proof is by induction on m; the lemma may be verified by explicit computation 
for fi and fg. Now assume the inductive hypothesis holds for f_m, . . ., fi, and use the 
equation (I2.2.3P to compute r_m-i term by term. Multiplication by E increases both the 
total and normal weight by 1 while leaving the ^-parity unchanged, and differentiation in 
x" does the same. Therefore \2C,\~^{—Er-m + dx^^m) has the properties we claim. The only 
remaining term in the expression for f_m-i is 

^ ^ ^ df{r,)D^{h). (3.1.3) 



2|e| ^ K\^ ^ ' 

-m<j<l 
-m<k<l 

\K\=j+k+m 
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The total weight of each term is (— j + 1) + {—k + 1) + K = m + 2, as required. Moreover, 
by the inductive hypothesis, the sum of the total weights, normal parities, and ^-parities of 
each term in fj and f^ is equal to zero mod 2 for all j and k. Therefore, the sum of the total 
weights, normal parities, and ^-parities of each term in ( I3.1.3p is i^ + i^ = mod 2. This 
completes the proof of the lemma. D 

The next step is to pass to the parametrix S{\). In this analysis, expressions of the form 

(iei-A)-'ierV(x,o (3.1.4) 

arise frequently; we define weights and parities of such expressions analogously to those for 
expressions of the form fl3.1.ip . Even though these expressions are restricted to x" = 0, they 
may involve normal derivatives of the metric, so normal weight and parity still make sense. 
A similar lemma holds for the Sj-. 

Lemma 3.1.5. For each m > —1, S-m{x, ^, A) may be written as a sum of terms of the form 
i\3.1.1\) : each summand has total weight m — 1. Moreover, the total parity, normal parity, 
and ^-parity of each summand are all well-defined, and their sum is equal to zero mod 2. 

Proof. The proof is extremely similar to the previous lemma, and again proceeds by induction 
on m. It is easy to check s_i explicitly. Now assume the inductive hypothesis. Using the 
previous lemma and the inductive hypothesis, the total weight of each term in the sum 
f l2.3.2p for s_m_i is {—k + 1) + (— j — 1) + /C = m, as required. Moreover, the sum of the 
total weights, normal parities, and ^-parities of each term in r^ and Sj is equal to zero mod 
2. As before, we conclude that the sum of the total weights, normal parities, and ^-parities 
of each term in fl2.3.2p is K + K = {] mod 2. This completes the proof. D 

We now pass from the parametrix S{\) to the heat invariants afc(x) themselves by using 
formula fl2.3.5p : fixing x, integrate in A and in ^. Each individual term in s_i_fc has the form 
f l3.1.4p . The contour integral in fl2.3.5p may be explicitly computed using Lemma [4.1.3j the 
result for each term is e~'^' times an expression of the form f l3.1.ip which has the same total 
weight, normal parity and ^-parity as the terms in s_i_fc. Now perform the ^ integral; this 
integral will vanish identically for all terms with ^-parity 1, and so each term which survives 
is a polynomial in the metric, its inverse, and their derivatives at x, with total weight and 
normal parity summing to zero mod 2. Since the total weight is A;, each term in ak{x) has 
normal parity equal to k mod 2. This completes the proof of Theorem 11.3.11 

3.2. Generalizations and applications. Recall that for the Laplacian on a closed Rie- 
mannian manifold, the local heat invariants are polynomials in the curvature tensor and its 
covariant derivatives |Gij . We may re- write Theorem 11.3.11 in this form as well. Let R^ be 
the curvature tensor for Vt and let H be the second fundamental form of M C fi; also let 
Vn and Vm be the covariant derivatives on Vt and M respectively. Then 

Theorem 3.2.1. For 1 < k < n — 1, the local Steklov heat invariant ak{x) may he written 
as a universal polynomial in the entries of the tensors V^Rn, < j < A; — 2, and V\.jH , 
0<j<k-l. 



Proof. The statement is an immediate consequence of Theorem 1.1.3 in |Gij . The restrictions 
on j follow from the weights in Theorem 11.3.11 D 

An interesting special case is when the embedding M G fl is product type: suppose that 
there exists a. 6 > and a tubular neighborhood ^ of M in O such that U is isometric to 
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M X [—5, (5]a;7i with the product metric. In this situation, the second fundamental form H is 
identically zero, as are all derivatives of the metric and its inverse in the x" direction. We 
claim 

Corollary 3.2.2. If the embedding M G Q is product type, the heat invariants ak{x) vanish 
for any odd k between 1 and n — 1. On the other hand, for even k between 1 and n — 1, the 
heat invariants are polynomials in the entries of the boundary curvature tensor Rm and its 
covariant derivatives V^jRm, < j < /c — 2. 

Proof. The vanishing of the odd heat invariants follows immediately from Theorem I1.3.1J 
each term in the universal polynomial for an odd heat invariant contains an odd number 
of normal derivatives of the metric, and hence contains at least one such normal derivative. 
Since the embedding is product type, each such term must be identically zero, and therefore 
the whole invariant is zero. 

As for the even heat invariants. Theorem 13.2.11 indicates that they are universal polyno- 
mials in Rq and its f2-covariant derivatives. Since all normal derivatives vanish, we may 
rewrite all of the f2-covariant derivatives of Rq_ in terms of M-covariant derivatives of Rm- 
Finally, examining the weights in Theorem 11.3. ![ the corollary follows immediately. D 

Note that in the product-type case, it is known that the Dirichlet-to-Neumann operator 
is a square root of the boundary Laplacian, in the sense that its square is equal, modulo 
infinitely smoothing operators, to the boundary Laplacian |Le] . 



4. Computation of heat invariants 

In this section, we compute and analyze the first few heat invariants, proving Theorem 
11.4.11 Throughout, we fix a point P on M and use the P-centered coordinates from Section 
2.1 to compute ak{P)- Recall from Section 2.1 that we may choose an initial orthonormal 
frame on M which diagonalizes the second fundamental form at P; we use this frame in all 
calculations that follow. 

4.1. Computations of oq and ai. We now evaluate the expressions (12.2. ip . f l2.2.2p . (12.3. ip . 
and f l2.3.3p at P using our local coordinates. These expressions, written out, are 



ri 



q-, = -^g-P^-^P; 



To 



^(E -iSi^tA.m) - ^ - ^a..* + f 



= i4 EE 9'"'(''"'^'' + ''■"'f °))(E '/'■'•ce) - A E 'r'-'a" 

7 a,l3 a,l3 ^^ „^^ 

+ ^ E^^^^''^^')^-" l°g^(^') + d.'^a^^i^'))^^ - \Y.3''^9.p,n. (4.1.1) 

Since we are using Riemannian normal coordinates at P, any first derivative in x" of the 
metric or its inverse, for 1 < a < n — 1, vanishes at P. Using f l4.1.ip and the fact that R is 
— 1 times the restriction of R to x^ = 0, we obtain that ri(P, ^) = |^| and 



2|^|2Z^ «- ^ 2 



a 
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Note that the subprincipal symbol of P, given in fl4.1.2p . has previously been computed by 
Taylor (Ch. 12, Proposition CI in |Tay] ) . We may also simplify (12.3. ip and f l2.3.3p : note 
in particular that any tangential first derivative of ri vanishes at P. Also remember that 
^°/3(p) = S'^^, gap,niP) = -2A„5„;3, aud ^"/'-'^(P) = 2XJ^i3. The results are 

Now plug these expressions into (12.3.50 . The contour integrals may be computed explicitly: 
Lemma 4.1.3. For any k > 1 and any a G M^, 

dX = —2Txi- 



lc{a-\f {k-l)V 

Proof. The proof is a simple computation with calculus of residues. D 

Using the contour integrals, then using the integrals in section 16.31 to evaluate the ^ 
integrals, we have, where Vn = Vol(S"~^), 

ao(P) = (27r)-"+i / e-l«l d^ = (27r)-"+V„r(n - 1); (4.1.4) 

a,{P) = i2n) l^^ I e '^'(1 ' ^^) ^^ = ,^,_iV9^^n-i 2^^° 



a 



2 7k-i ' \i? ' 2(n-l)(27r) 

Vn{n - 2)V{n - 1) 



, , , Hi. (4.1.5) 

2(27r)'"-i ^ ^ ' 

This completes the proof of the expressions for oq and ai in Theorem 11.4. 1[ 

4.2. Computation of 02- The computation of 02 is somewhat more involved; this is typical 
of heat invariant calculations, which tend to increase dramatically in complexity as one goes 
farther out in the expansion. As a starting point, we write out the expression for s_3 given 
in (ESS): 

S--AP, e, A) = {r,{P)fm - X)-' - (r_i(P))(|e| - A)-^ - .(%i ■ d^voMl - Xy' 
+ (%! ■ d^d^r, . d^r,)m\ - X)-' + (n - X)-\ J2 ^{dfr,){P){d^s_,){P)). 

\K\=2 

Simplifying, using the fact that first derivatives in x of ri, and therefore also first derivatives 
in X of s_i, vanish at P: 

s,s{P, e, A) = (ro(P))2(|e| - A)-3 - (r_i(P))(|e| - A)-^ - ^i^^r, ■ d,r,m\ - X)-' 
+ (5^(9eri9^.n9..9..r0)(|el - A)-^ - ( J2 ^id^n)iPMr,)iP)m\ - X)-'. 

7,e \K\=2 

As before, plug this expression into (I2.3.5P and evaluate the contour integrals. We also switch 
from rj to fj for later ease of computation. After all this, we obtain the following expression 
for 02: 



a2(P) = (27r)-"+i / 6(P, Oe-'?' d^ (4.2.1) 
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where 

- I E(^en%ri9..9..r~0 -^{J2 ^id^h)iP,md^h)iP,0)- (4.2.2) 

7.f \K\=2 

Our strategy is a direct approach: write out b{P, ^) in terms of the metric and then 
integrate to get 02 (P). In this calculation, b{P,^) splits naturally into two components 
bn{P,C,) and bt{P,C,), which we call the normal and tangential components respectively. To 
define 6„ and bt, first write out the expression for f_i in terms of fi and fo: 



'^' ^ - |i^|=2 



iro)^ - I 



+5x"ro + ^ J] ^"''(^a/^.nro ] . (4.2.3) 



The normal component 6„ consists of the terms involving normal derivatives of the metric; 
as we will see, these are the first term in f l4.2.2p along with the restrictions to x" = of the 
first term and the last two terms in (14.2.31) . The tangential component bt is the remainder. 
The expressions for 6„ and bt may be written out and simplified; note in particular that 
the term in f_i involving a first derivative of the metric, dxi'Ti, is zero at P. Writing out 
the multi-index notation, and noting that the factor of K\~^ precisely compensates for the 
double-counting, we have: 

bniP,0 = il + ^)^~o(^'0 + ^(5x."ro)(P,0 - ^E^"^~o(P,0; (4.2.4) 



bt{P,0 



-i(,7^ + ^777) ^ d^jTi ■ (9^7 ro - ^ X] ^^"'^^ ' ^^''^^ ' ^^''^^''^^ 



2 ' 2\iy^^'^ ^'^ 6 

' ' 7 7,<: 



(P,0. (4.2.5) 



~^4 + 4^^ ^ d^.d^.h ■ d,-,d,.h 

Now we write bn and bt in terms of the metric. Recall that when integrating over the 
tangent space, any term which is odd in ^, or even odd in ^° for any particular a, will 
integrate to zero. We therefore define an equivalence relation = on functions of ^ by writing 
a(0 = &(0 iff 






Lemma 4.2.6. Let Rij be the Ricci tensor on M. Then: 






^|^|3^^ v^.'2|e|^"" 8|e|^^""""' 41^1^-— I^P 
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m^,o = (y^ + y^)E^-(^ 



a\2 



The proof is a direct calculation; we plug in the formulas for the fj and simplify. The 
analysis of bt uses the Taylor expansion of a Riemannian metric in normal coordinates as 
well as symmetries of the Riemann curvature tensor. The details of the proof are deferred 
to section 16.11 

Since we have found 6 = 6„ + 6f , we proceed to compute 02 (-P), which we correspondingly 
write a2,n{P) + 02,t(-P)- In the computation, the integrals in section [673] are useful. Each 
integral contains a factor of V"„, and since each integral has either k = or k = —1, we 
can also bring out a factor of r(n — 2) (and then multiply the k = integrals by n — 2 to 
compensate). We obtain: 

V^Tjn - 2) ^ 

«2,i(^) = -.^/^ x^_i Rm, 



12(27r; 



n-l 



a2,n{P) 



VnTjn - 2) 
(27r)"-i 



iEA«Mn-l)(l-;^ + i±^) + ii:A.A 



1 + 26^, 



- a, 13 



n^ 



'a a a a,l3 



8(n-l) 
Simplifying and using f l2.1.ip . we find 



(4.2.7) 



a2,n{P) 



VnT{n - 2) 



K27r; 



71-1 



/ r- 2 n + 3 , >r^ , , 

(n - 5 + 7 + — j) 2^ A„A^ 

a,/3 



n 



+ ( 



1 

,2n + 6 
n^ — 1 



n^ 



+ 4 



n 



1 ^-^ n — 1 ^-^ 



Finally, recall that 



[n 



l?Hl- Y.^l = {n- If Hi - (n - l)(n - 2)if2, 



(4.2.8) 



(4.2.9) 



",/3 



so we may rewrite the eigenvalue sums as combinations of iff and if 2- Moreover, the last 
term in the expression for a2,n may also be rewritten in terms of mean curvatures: 



Lemma 4.2.10. 



2_^9aa,nn{P) — Rf 



R 



M 



2{n - lyHf + 3{n - l){n - 2)H2. 



This proof is deferred to section [6721 From Lemma [4.2.101 and (I4.2.9p . we conclude after 
some algebraic manipulations that 



a2,n(P) 



VnT{n - 2) 



K2vr; 



n-l 



(^_l)(^_2)(n2-n-4) 2 
^1 



n 



n+1 
Combining this with 02,4 (P) yields Theorem 11.4.11 



n{n-3){n~2) n-2 ^ 

--"2 + Ti-^n — n-M) 



n 



. (4.2.11) 
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4.3. Example: heat invariants of balls in M.^ and W^. In this section we verify the 
formulas for the heat invariants obtained using Theorem 11.4.11 by calculating them directly 
for the unit balls B^ C M^ and B^ C M"*. It is well-known that the eigenvalues of the Steklov 
problem on a unit ball B" C M" are given by a sequence of natural numbers k = 0,1,2,3 ... , 
with each eigenvalue k repeated according to its multiplicity 

{2k + n-2){k + n-3)\ 



rrik 



kUn-2) 



^^-°-E '""".;"^""^" ^- (^-3.1) 



Note that the eigenf unctions of the Dirichlet-to-Neumann operator V on S""^ are spherical 
harmonics, and the numbers ruk are the multiplicities of the Laplace-Beltrami eigenvalues 
on S"~^ (see, for example, |Po] ) . The Steklov heat trace on a ball is given by an explicit 
formula 

{2k + n-2){k + n-3)\ 
k\(n-2)\ 

Let us compute its asymptotics as t — )■ 0+ for n = 3 and n = 4. For n = 3, the series (14.3. ip 
takes the form 

£(2* + 1)'-" = 4^2 = '"' (2 + « + T + 0(1")) , 

k=0 \ J \ / 

and therefore the corresponding heat invariants are oq = 2, ai = 1 and 02 = 1/3- Taking 
into account that V3 = 27r, Hi = H2 = 1 ai each point x G S^ and Vol(§^) = Ait, we obtain 
precisely the same values for oq, ai and 02 from Theorem II. 4. 1[ 
Let now n = 4. Then (14.3. ip takes the form 

CO 00 2tf-\ j_ t\ 

J2ik + Ife-"' = e* ^(A; + 1)^6'^'+'^' = "^ j _ ^^ = r^ (2 + 2t + t^ + 0(1')) . 

k=0 k=0 ^ ' 

The corresponding heat invariants are a^ = 2, ai = 2, 03 = 1. These results are in agreement 
with Theorem II. 4. H as one can easily check, taking into account that V4 = 47r, Hi = H2 = 1 
for any x E S^ and Vo\{S^) = 2%'^. 

5. Proof of Theorem 11.5.31 

In this section, we adapt an argument of Zelditch in fZ] and combine it with an analysis 
of the first three heat invariants to prove Theorem 11.5.31 

5.1. Zelditch's theorem on multiplicities. In order to make use of the heat invariants, 
it is first necessary to prove that the boundary of any domain in M^ with the same Steklov 
spectrum as a ball must be simply connected. This follows from a stronger result, which is 
an adaptation of a similar result of Zelditch for the Laplacian (Theorem A in [Z]). 

Proposition 5.1.1. Let U be a Riemannian manifold of dimension n with smooth boundary 
{Y,g), and let = Xq < \i < X2 < . . . be the distinct eigenvalues of the Dirichlet-to- 
Neumann operator V, with multiplicities rriQ, rrii, .... Suppose there exists a > such that 

mk = ak''-^ + 0{k''-^). 

Then {Y,g) is Zoll; that is, the geodesic flow, which we call G*, is periodic. 
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Proof. The proof is closely analogous to the proof in [Z] , which is in turn based on work of 
Ivrii and Hormander. Let Il*{y,r]) be the niicrolocal period function on T*Y \ {0}, equal to 
the period of the geodesic corresponding to {y, rf) if it is periodic and equal to oo if it is not. 
Assume for contradiction that there is a non-periodic geodesic; then the set Vt of (y, rf) for 
which n*(|/, 77) > T is a nonempty open cone for each positive T. Fix a large positive T to be 
determined later. As in [Z], let S be a self-adjoint pseudodifferential operator of order zero 
with microlocal support nontrivial and contained in Vt-, and let h be the principal symbol of 
B*B. Then define 

N{X,B*B) = J2 *TrB*B\E„ 

Afc<A 

where i?fc is the P-eigenspace corresponding to A^, and the sum ^* is over distinct eigen- 
values. 

There is again a modified Weyl law for N{X,B*B) given in the proof of Theorem 29.1.5 
in Hormander |Hoj . The difference with [Z] is that the subprincipal symbol of P, rQ{y,rj), is 
not necessarily zero. On the other hand, the subprincipal symbol of B*B is still zero, and 
the integral of the Poisson bracket {b, \ri\} is still zero, so from [Ho] we have: 

N{X, B*B) = (27r)-("-i) ( ff hdydr, + d^ If r,{y, r])b{y, r]) dy dr]] + R{X, B*B), 

\JJ\ri\<X JJ\ri\<X / 

(5.1.2) 
where the remainder -R(A, B*B) satisfies the estimate 



limsupA-("-2) \R{X,B*B)\ < ^ 



h dy drj 

\v\<i 



(5.1.3) 



'^'-^-^K- 



N{X, B*B) < ^^"^J" ^h J2 *^r' + g'T(l). (5.1.4) 



We write b = Jj ^b dy drj and note that since B has nontrivial microlocal support, 6 > 0. 
The first two terms in (I5.1.2P are continuous in A, which allows us to write 

TrB*B\E, = \im{R{Xk + e, B*B) - R{Xk - e, B*B)). 
Hence for A^ greater than some Xq{T) (it may depend on T), we have 

TrB*B\E, < '^^'^^'^h xi 
Summing gives 

T 

Afc<A 

Next, use the multiplicity assumption, precisely as in [Zj, to rewrite (I5.1.4p : the proof is 
identical to that in [Z], albeit with A^ here replacing -y/Afc in [Z], so we omit it. We conclude 

that 

iV(A, B*B) < ^bX"-^ + C(A'^-2) ^ 0^(1)^ (5_X_5) 

where Z) is a constant depending only on n and the volume of Y . Note that the first 
term in (15.1.2^ is exactly (27r)^*^"^^^6A"~^, and all other terms are (9(A"~^). Thus fixing 
T > (27r)"-iD and then letting A -^ 00 in flET3D contradicts flHX^ . 

We have now shown that all geodesies on Y are periodic; as in [Zj, we now apply theorem 
(0.40) in |Be] to conclude that all the closed geodesies must have a common period, and 
therefore that (F, g) is ZoU. D 
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Remark 5.1.6. Suppose, on the other hand, that f/ is a manifold with boundary Y having 
the property that the periodic geodesies on Y form a set of measure zero in T*Y. Then, as 
follows from the results of jDG] (see also [Hoi Chapter 29]), there is a two-term Weyl law 
for the Steklov eigenvalues. We may write it as 

where a^ and ai are the heat invariants. Here we have used the Laplace transform to relate 
the heat trace and the counting function. 

Now assume that that fi is a compact domain in M^, with smooth and connected boundary 
M, which has Steklov spectrum equal to that of the ball of radius p in M^; this spectrum is 
given up to scaling in section 1^^ and satisfies the hypothesis of Proposition 15. 1 . 11 Therefore 
any domain in R" which is Steklov-isospectral to the unit ball must have Zoll boundary. 
Since all connected Zoll surfaces embedded in M^ are topological spheres [Be], we conclude 
that the boundary M is simply connected. On the other hand, since there are large families 
of Zoll surfaces |Guj . we need more information to conclude that M is a sphere. 

5.2. Application of heat invariants. We now use the heat invariants we have computed 
to finish the proof of Theorem 11.5.31 By the previous paragraph, the Gauss-Bonnet theorem, 
and f ll.4.8p . we know that the second Steklov heat invariant 



^2 = 

IGvr 



L//?-^>')-ii//?-^^- (^■^■^; 



Therefore j]yjHf = j' Hf. On the other hand, we already know from the first two heat 
invariants that Vol(M) and /^j-f^i are Steklov spectral invariants, so Vol(M) = Vol(S'p) and 
/m^i = /5„^i- Therefore 



^Yo\{M)(j Hi 



1/2 
2 



Hi 

M 




Hi 

So 



0. (5.2.2) 



By the Cauchy-Schwarz inequality. Hi must be constant on M. However, the only compact 
surfaces of constant mean curvature in R^ are spheres [Alj, so we conclude that M is itself 
a sphere of radius p and therefore Vt is isometric to Bp. This completes the proof. 

6. Proofs of auxiliary lemmas 
The following two proofs are given by lengthy but straightforward computations. 

6.1. Proof of Lemma 14.2.61 To analyze the normal integrand (14.2.41) . first recall the ex- 
pression for fo(P, ^) given in (14.1.21) . Plugging in this expression takes care of the first and 
last terms in (I4.2.4P : it remains only to analyze the middle term i^^dx^rQ. Examining the 
expression (I4.1.ip for fo, we notice that the first and third terms (the imaginary part of fo) 
are odd in ^, and hence the same is true after applying a normal derivative and dividing by 
2|^|. Therefore 



wf-'''^h'A-w^/"''°''-\^f ''-■'■' 
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^ ^aa,nn^^a)2 + Y^ ^2 _ ^ 



^i^i3z^- - . 2iei^ " m 



Qaa, 



nn- 



In the final step of tliis calculation, we used the fact that |^|^^'C"'C^ — whenever a 7^ /3. 
Combining this calculation with the rest of the integrand completes the proof of the normal 
portion of Lemma 14.2.61 

Tangential integrand: initial computations. Now consider the tangential portion. The first 
component of 6t is — 2(| + 2^1) Z^^ d^-yrid^i'To. Notice that 

which is odd in ^. Each term in the real part of d^-i'To will be a power of |^| times a polynomial 
of even degree in ^, so after multiplying by ^'''/|^|, it will be equivalent to zero. We therefore 
need to consider only the imaginary part of dxir^, and hence only the imaginary part of fo 
itself. However, from ( 14.1.ip . each term in the imaginary part of fo is multiplied by a first 
derivative of the metric or the log of the volume element, which is zero at P. Therefore, 
the dx-i derivative must always hit that term, or else the result is zero at P. Relabeling the 
dummy variable 7 in ( 14.1. ip as e, to avoid confusion with d^i, we conclude that 



a.-oW=(-^E 



. 01,13 a, (3 



J2(l9'''^K-\og6 + g-^''^^)A (P). (6.1.1: 



21^1 ^'2^ d^.dx 

1^1 a,f3 ^ ^ 

Now plug in the values of the metric at P, and then multiply by i^'^{\ + ^X As before, we 

also note that any term of the form |^|~'^^"^^ is equal to zero whenever a 7^ /3. We get that 
the first component of bt{P,^) is equal to: 






a\2 



For the remaining terms, we first compute second derivatives of f 1 in both the x variables 
and the ^ variables: 

(a,A-r.W) = (^^p ^(P) = j^-Yi. 
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From fl4.2.5p . the remaining terms in bt{P,^) are hence equal to 

(y^ + ^ + ^) Y:^3-'-^{p)cecc - (^ + ^) j:^9-''''{p)ce- (6.1.3) 

Curvatures. We now use the well-known Taylor expansion of the metric in Riemannian 
normal coordinates (see |Vi] , for example) to relate the second derivatives of the metric and 
the volume element to the intrinsic curvatures of the boundary M. Let Rijki be the Riemann 
curvature tensor of the boundary M at the point P; then we have 

An easy inverse argument using Taylor series gives: 

g'^'^ix) = 6^p + ^-R^^p.x^x'' + 0(|xn, 

and hence 

^a/3,7.(p) = 1 (p^^^^ + R^^^^y (6.1.4) 

As for the volume element, we have from fVil: 



5{x) = 1-\r^,x^'x'' + 0{\x\^), 

so by elementary Taylor series arguments 

log(5(x) = -^i?^,x'^x^ + 0(|xp), 

and therefore 

(%aelog5)(P) = -^{R,, + R,^) = -'^R,,. (6.1.5) 

Combining these observations with f l6.1.2p and f l6.1.3p . and also removing the a ^ [5 part 
of the last term of (16.1.30 (which is equivalent to zero), we get: 

111 



1 



^)5^i?a,.,(r)'. (6.1.6) 



Final simplification. To simplify (I6.1.6P further, we use the symmetries of the curvature 
tensor: 

^ijkl ^jikl ^ijlk ^klij- 

From these symmetries, we immediately conclude that Raapi^ = for all a and /3, and also 
that Rajspa = —Rai3ai3- Heuce the third and fourth terms of (I6.1.6P cancel. Moreover: 

Proposition 6.1.7. For any k > —3 — n (so that the integral over the tangent space makes 
sense), 
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Proof. Integration of a multiple of ^"^^^'^^'^ against an even function of ^ vanishes unless 
the four indices pair off into groups of two. On the other hand, if all four are the same, the 
curvature coefficients are zero. We may therefore consider each possible pairing separately 
without worrying about double counting; for example, ii a = /3 and 7 = e, the sum becomes 

Relabeling and adding up the possibilities, using the symmetries of the curvature tensor, the 
first sum becomes 

a,l3 

which completes the proof of the proposition. D 

As a consequence of the proposition and the preceding remarks, the first term in (16.1. 6p 
vanishes, and we are left with 

This completes the proof of Lemma 14.2.61 

6.2. Proof of Lemma 14.2.101 The proof is by computation of the scalar curvature R^. At 
the point P, using Einstein notation, we have 

Rn = '^aa,c ~ '^ac,a + ^ acX cd ~ ^ aJ^ ad- (6.2.1) 

Moreover, if the sums are taken only over the indices from 1 to n — 1, rather than from 1 
to n, we obtain the scalar curvature of the boundary. So Rn equals Rm plus the terms in 
(16.2. ip where at least one of a, c, or d is n. On the other hand, the Christoffel symbols are 
given by 

m 

Recall that Qan and g""^ are identically zero in a neighborhood oi P dVt for any a < n — 1, 
and gnn and (/"'" are identically 1. As a consequence, whenever two or more of i, j, and k 
are equal to n, the Christoffel symbol Ff • is zero in a neighborhood of P. 

We first analyze the first two terms in (16.2. ip where at least one index is n. Using the 
observation above on the vanishing of the Christoffel symbols, we obtain that at P, 

re— 1 n— 1 n— 1 n— 1 

E-pn \ "^ -pc \ ^ -pn \ ^ -pa 
aa,n / j nc,n / j aa,n / j nct,n 

a=l c=l a=l a=l 

^ 71—1 _ n—l 

7) / J Or" \ / J \9am,a ~r (jma,a 9aa,m)9 ) ' T\ / j ^x^-y / ^ \9am,n ~r (Jmn,a 9na,n)9 ) 
a=l m 7=1 "I 

^ n—l ^ n—l n—l n—l n—l 

^ 7j / _, [ 9aa,nn) ^ / _, '^x" ( / ^ 9am,n9 ) ^ / ^ 9aa,nn + -^ / ^ ^q.- 
a=l a=l m=l a=l a=l 

Now consider the final two terms of (I6.2.ip . Since there are no derivatives of the Christoffel 
symbols involved in these terms, we can plug in the metric at P and write 

r-,(^) = ^{9jk,i + 9ki,j - 9^J,k){P). 
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If more than one of i, j, or k is n, then F^- vanishes; on the other hand, for a,P ^ n, we 
compute 

If more than one of a, c, or d is equal to n, the final two terms of fl6.2.ip vanish; if none are 
n, then we get part of the scalar curvature Rm- If one is n, then there are three possibilities: 

• Suppose a = n, c ^ n, d ^ n. Then the third term of (16.2.11) is zero; the last term is 

71—1 n—1 



$:-MA.f=-E^- 



c,d=l a=l 

• Suppose a y^ n, c = n, d j^ n; the third term is again zero, and the last term is 

n— 1 n—1 

a,d=l a=l 

• Finally, suppose a j^ n, c j^ n, d = n. Then we get: 

n-l 

J2 Aa(-Ae) + Y. -i-^aSac)iXaSac) = "(^ " ifH^, + ^ A^ . 
a,c=l a,c a 

Combining these three computations with the first two terms of (16.2. ip yields 

Rn = Rm + 3 2_^ Aq, — 2_^ 9aa,nn " (^ " 1) H^. 
a a 

Rearranging and using (14.2.90 completes the proof. 

6.3. Auxiliary Integrals. 

Lemma 6.3.1. Let a and /3 be any two distinct integers between 1 and n—1. Moreover, let 
Vn = Vol{S^^^). For any real k such that the following integrals converge, we have 

e"'^'|eNe = Kr(A; + n-l); (6.3.2) 



f e-l«l |e|'~'(r)' d^ = -^r(A; + n - 1); (6.3.3) 

Jk„-i n-l 

[ e-^e-\CM = -^nk + n - 1); (6.3.4) 

J^n-i n^ -1 

[ e-\^\\^\'-\e)\ffd^ = -^^{k + n-l). (6.3.5) 

J^n-i n^ - 1 

The proof is a direct calculation in n-dimensional spherical coordinates; the resulting 
trigonometric integrals may be evaluated using, e.g., |GR] . 
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